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We explore the phase diagram of NJL-type models near the chiral critical point allowing for 
phases with spatially inhomogcneous chiral condensates. In the chiral limit it turns out that 
the region in the mean-field phase diagram where those phases are energetically preferred 
very generically reaches out to the chiral critical point. The preferred inhomogeneous ground 
state in this vicinity possibly resembles a lattice of domain wall solitons. This raises the 
question of their relevance for the phase diagram of QCD. 

I. INTRODUCTION 

The phase diagram of quantum chromo dynamics (QCD) as a function of temperature T and quark 
chemical potential /x is expected to exhibit a rich phase structure In particular the nature of the 
chiral phase transition, the location of the chiral critical point and the properties in its vicinity have 
attracted a lot of interest. Experimentally, this exploration is one of the main research goals at the 
future FAIR facility in Darmstadt 0] • Theoretically, it is heavily investigated in phenomenological 
models whereas ab initio calculations are still limited to small values in /i/T. Depending on the 
approach and/or details of the model a great range for the possible location of the critical point 
has been found 

Generic, model-independent scenarios near the critical point can however in principle be dis- 
cussed via a Ginzburg-Landau (GL) expansion. The idea is to expand the effective action of the 
order parameter and to obtain an expression for the thermodynamic potential as a functional in the 
order parameter. For inhomogeneous phases, characterized by a spatially varying order parameter, 
this can be generalized by a combined gradient expansion. Treating gradients and the magnitude 
of the order parameter to be of the same order, a generic expression for the GL functional in the 
vicinity of a second order phase transition for a theory with real order parameter 0(x) invariant 
under (/'(x) — > — (/'(x) takes the form 



r?GL(T, ^; </.(x)) = C2(T, p)(t){^f + C4,a{T, li)4>{x)^ + Ci,b{T, /u)(V0(x))2 

+C6,a(T, ^)0(x)6 + Ce,b{T, ^)(V0(x))20(x)2 + C6,c(T, /i)(A,/.(x))2 . (1) 

Considering a periodic ground state with a Wigner-Seitz cell V, the thermodynamic potential 
difference between symmetric and spontaneously broken phase is then given by AO = y Jy Qqi . 
Furthermore we have to require cg^j > in order to have a bounded potential. For homogeneous 
phases with 0(x) = (pQ a second order phase transition is then occuring at C2 = 0, C/^^a > 0, 
which turns into a weak first order phase transition at the critical point where c^^a changes sign. 
In addition inhomogeneous phases can become energetically favored when C4^b < 0, since small 
curvatures then lead to a gain in free energy. 
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The main idea of this work is to perform such an expansion within an NJL-type model on 
mean-field. For simplicity we limit ourself to the chiral limit and do not include [/A(l)-breaking 
terms. Furthermore constraining ourself to condensates of the form {tpi{x)tl^j{x)) oc (f){x)6ij, where 
i,j label flavors and chiralities, the residual global symmetry is Z2 and the GL functional should 
take the form as stated in Eq. ([1]) . Assuming a sensible regularization procedure we find c^^a = C4,b 
and related to that an inhomogeneous phase starting at the critical point. 

This finding is actually very similar as in the one-dimensional^ Gross-Neveu (GN) model, which 



is renormalizable and where the mean-field problem can be solved analytically 0, l5[. In the mean- 
field phase diagram of the GN model a symmetric, a homogeneous dynamically broken and a 
inhomogeneous dynamically broken phase exist and meet at the chiral critical point. A similar 
behavior persists also for finite quark masses 

To conclude the introduction we want to clarify, that we do not necessarily want to conjecture 
precisely this scenario for the phase diagram of QCD. Rather we would like to suggest that in- 
homogeneous phases may play a more important role than anticipated. It should be noted that 
the NJL model shares the global symmetries with QCD and is therefore a reasonable candidate 
to study the phase properties near the critical point even if the location of the critical point may 
differ. So far these chiral crystalline phases have mainly been discussed at vanishing tempera- 
tures 0, S, S)E3l- Here the discussed scenarios are qualitatively different in that the modulations, 
i.e. the magnitude of the gradients, are not small but typically of the order of the Fermi momen- 
tum. The relevant dynamics of particle-hole scattering for one plane wave is then concentrated on 
a resticted vicinity of two antipodal points on the Fermi surface. An exception to this is the dual 



chiral-density wave which is sufficiently simple to solve the mean-fie ld p roblem On a tech- 

nical level it has strong similarities with the Fulde-Ferrell phase Il3| discussed in the context 
of (color-)superconductors. The latter is known to be technically simple but disfavored compared 
to other inhomogeneous phases Q, [3]. In particular the obtained order of phase transitions 
when going from homogeneous dynamically broken to inhomogeneous to unbroken phase might be 
misleading [l^ . 

II. THE GENERALIZED GINZBURG-LANDAU EXPANSION 

We consider NJL-type models of the form 

C = ^i-f^d^^P + Gii^ipf + . . . , (2) 

where ip is the 4A^jA^c-dimensional quark spinor for Nj flavors and Nc colors, are Dirac ma- 
trices and G is called the scalar coupling. The dots indicate terms that vanish in our mean-field 
approximation. Allowing a mean-field value (VS(x)^/>(x)) = — 2^M(x) for the chiral condensate, 
the mean-field Lagrangian density takes the form 

Cmf = i^ii^d^ - M(x))^ - . (3) 



^ In this work we only refer to spatial dimensions. 



3 



In the case of a periodic condensate with Wigner-Seitz cell V and using the imaginary-time for- 
malism (see e.g. Refs. [1]), we therefore obtain for the mean-field thermodynamic potential 



n{T,fx) = -^In / V^Vi^expl [ (^AfF + /^V^tV) 1 

^ J \Jxelo,^]xV J 

= -yTTn,cj,v Log {S-') + - + const. , (4) 



with inverse propagator 

S-\x,y) = (i7^a^-Af(x))5(4)(x-y) (5) 

and the functional trace acting on the direct product of Dirac, color, flavor and coordinate space. 
Since the functional logarithm can only be evaluated in special cases such as the dual chiral- 
density wave we now expand in M(x). Substracting the leading order corresponding to the 
thermodynamic potential of the unbroken phase, we formally arrive at 

/^) = E ^TrD,c,/,y {SoMT + y [ ^ • (6) 

n>0 

Here we have introduced the bare propagator 5*0 = S\]^^(^^■^^o and a short hand notation for 
TrD,cjy{SoM)'' = [ f ...[ Tr^.c,/ (M(x)5o(x, a;2)M(x2) . . . M(x„)5o(x„, x)) . 

J X J X2 J Xn 

(7) 



The domain of integration for x is [0, X y and [0, ^] x for X2, ■ ■ ■ ,Xn- In the chiral limit the 
expressions for odd values of n vanish. Furthermore we can expand the condensate around x as^ 

M(xO = V i-D"M(x)(x„-xr (8) 
^ — ' a! 

|o|>0 

and can extract the GL functional to any desired order in gradients and order parameter. Neglect- 
ing possible issues with the regularization for the moment, we can go to momentum space using 
So{x,y) = TJ2n (2tt)'-^ {PiJ.l^)~^ ex.'p{ip{x — y)) , where pQ = Un = (2n -|- l)7rT. It is then a tedious 
but straightforward exercise to work out the explicit expression stated in Eq.([T]). All coefficients 
result in the evaluation of similar integrals, only 

takes a slightly more complicated form, in that the integrand involves total derivatives. In addition 
part of the expressions are of course formally divergent. 



^ We are using multi-index notation. 
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Since our model is non-renormalizable the divergences cannot be absorbed into a redefinition of 
tfie coupling. Instead usually a regularization, as part of the phenomenological model, is introduced. 
Due to this ad hoc procedure a generalization of the regularization to inhomogeneous phases is often 
not unique^ as already discussed in the context of inhomogeneous color-superconductors [3] • We 
could therefore take the pragmatic viewpoint that a generalization of any such ad hoc regularization 
procedure to inhomogeneous phases is assumed to be such that the total derivative terms vanish. 
An alternative approach is a regularization scheme that does not rely on an homogeneous ground 
state, e.g. a propertime regularization for the functional logarithm in Eq.(j3|). In this case it is 
possible to show that the total derivative terms strictly vanish. For both of these viewpoints we 
can define the regularized expressions 

(fp 1 , 82 



On 



{-DUn^n^tY^I tC^t-— -W^ + I^ (10) 



and (not replacing M(x) by (/>(x) for simplicity) finally obtain 

J^GL(r,^;M(x)) = ^M(x)2 + ^ (M(x)4 + (VM(x))2) 



^ !V1 \ X I -t- I I \ / /I /I < 1 > 

.'^ ( M(x)6 + 5(VM(x))2m(x)2 + i(AM(x))2 
6 V 2 



(11) 



Interestingly the ratios of the prefactors within each order of the GL expansion turn out as those 
in the one-dimensional analogue 0, 0]. This is not the case in superconductors 17 1 where the 
underlying dynamics is different, namely coming from particle-particle and hole- hole scattering 
near the Fermi surface instead of particle-hole scattering in presented case. 

It is now in principle straightforward to choose a regularization scheme and adjust the model 
parameters to QCD phenomenology. Since it is however known that differing choices of regular- 
izations and model parameters lead to a spread in the location of the critical point (see e.g. the 
collection in Ref. 0]), we prefer to make a more qualitative statement and keep the discussion on 
the level of the GL coefficients Although the model may not have the same critical point 

as QCD, it should be in the same universality class and therefore show a similar behavior in the 
vicinity of the critical point. 



III. INHOMOGENEOUS GROUND STATES 



Given the GL functional in Eq. pij) and some critical point where a2 = 04 = and > 0, we 
can explore the phase structure in its vicinity. Limiting ourself first to homogeneous phases with 
M(x) = Mq, we have the following known behavior: 

2 

04 > 0: Second order phase transition at 02 = 0, where Mq ~ — ^ and AQ = —477- for 02 < 0. 



^ This applies in particular for the widely used three-momentum cutoffs/formfactors. 
For all conventional choices the presented model has a critical point in the phase diagram where 02 = Q4 = and 
ae > 0. 
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FIG. 1: Pictorial presentation of the phase diagram in terms of the GL coefficients: The gray domain 
corresponds to the homogeneous dynamically broken ground state, the shaded gray to the solitonic ground 
state (at least when restricting to one-dimensional modulations in the order parameter), whereas in the 
transparent domain the unbroken phase is preferred. A is an arbitrary scale. Also stated are various lines 
discussed in the text. 



04 < 0: First order phase transition at 04 = — ^^02^6 where the broken solution has Mq = 
~iof ■ The dynamically broken solution continues to exist as a local minimum up to a4 = 
— -v/4a2a6- 

Allowing for inhomogeneous phases we may expect a spatially varying order parameter for a4 < 
since here small curvatures can lead to a gain in free-energy. Even within the GL approximation the 
determination of the ground state is not straightforward as we need to minimize a non-quadratic 
functional. Focusing on one-dimensional inhomogeneities, i.e. M(x) = M(z), the solutions to 
■^AQ = are actually known from the investigation of one-dimensional models E^- They 
are expressed (up to an arbitrary shift) in terms of the elliptic Jacobi sn-function as 



(12) 
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where v G [0,1] and q being a scale related to the maximum of Mx^i^z) and the extension of 
a soliton in the chosen z direction (both scales are related in our case). For = 1 we have 
^\d{,z) = gtanh((/x), i.e. a single soliton and for — > the shape becomes more and more 
sinusoidal albeit the amplitude also goes to zero. From previous investigations 0] it is known that 
when increasing 02 from zero we reach a second order phase transition into an inhomogeneous 
phase with q = Mq and v = 1. At this point the free-energy of a single soliton becomes negative 
leading to its formation. By using Mq known from above and checking where ^Ar2|7y^(x)=Mio(z) 



changes sign at = 0, we obtain 04 = — y ^0205 for this point. We arrive at the onset of infinitely 
far separated solitons. Further increasing 02 decreases v until it reaches zero. Since q stays finite 
the overall magnitude of Mi£){z) given by y/vq then vanishes and we find a second order phase 
transition to the unbroken phase. 

In case of a second order phase transition from the inhomogeneous phase to the unbroken phase, 
the value of in terms of a2aQ is actually universal also for higher dimensional modulations of 
the order parameter. Since in this case M(x) is parametrically small, we can neglect non-quadratic 
terms in the GL functional. Consequently the variation ^^Afi leads to a linear partial differential. 
We can then optimize the value of 04 by varying the momentum q of the Fourier components of 
M(x) and find 04 = — Y^|a2«6 foi' the transition line where |q| = 

We do not want to address the general question whether an inhomogeneous phase with a higher 
dimensional modulation could become favored in the vicinity of the critical points, but it may very 
well be that the one-dimensional modulations are generally preferred there as numerically confirmed 
in Ref. [3] for the analogous case of inhomogeneous phases in paramagnetic superconductors. 

We summarize our findings in Fig. [1] by illustrating the ground states in the GL coefficients' 
phase diagram. Obviously the inhomogeneous phases modify the region where a first oder transition 
is expected when restricting to homogeneous ground states only and replaces this transition by two 
second order phase transitions. 



Starting with generic NJL-type models in mean-field approximation, we have preformed a gen- 
eralized GL analysis in the vicinity of the chiral critical point. Since the order parameter and 
the gradients of the modulation are parametrically small, the expansion is arbitrary close to the 
mean-field description. Requiring a sensible regularization scheme the results are astonishingly 
similar to results in the one-dimensional GN model. Consequently, we find two second order phase 
transitions from homogeneous to inhomogeneous dynamically broken to unbroken phase, instead 
of a first order transition from broken to unbroken phase. 

Various questions can still be addressed: It would be interesting to actually compute the GL 
coefficients in a regularization scheme with parameters adapted to QCD phenomenology and to 
estimate the extension of inhomogeneous phases in the phase diagram of NJL-type models. This 
should be combined by a complete mean-field study when at least allowing for a one-dimensional 
inhomogeneity in order to check the range for the applicability of the GL expansion. It should 




IV. DISCUSSION 
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also give a different phase diagram than obtained for the analogous one-dimensional models, where 
the inhomogeneous phases reach out to infinitely high chemical potentials. Another interesting 
question here is whether the order of the phase transitions change when going away from the 
critical point as recently found at vanishing temperatures for (color-) superconductors [3]. Also 
the competition with color-superconducting phases would be interesting as the latter, at least in 
such model studies, should not reach out to the critical point. 

There are also various ways the GL analysis can be improved: It would be interesting to include 
finite and mutually differing current quark masses and to check how the presented scenario is 
affected. This would in addition require the use more complex order parameters {tpitpj), where the 
indices i, j should include flavor and chirality. Also more complicated interactions or models could 
be considered, such as vector interactions, the t'Hooft interaction when in particular including 
the strange quark or the Polyakov-loop NJL model. Incorporating time derivatives one may try 
to explore dynamical properties of fluctuations and discuss possible consequences for heavy ion 
collisions. 

More complicated but very interesting problems would be whether those phases could be ex- 
plored beyond mean-field approximation, how the relation between the GL coefficients change when 
including higher fluctuations and whether the phase diagram of QCD may include those or similar 
phases, in particular in the region close to the chiral critical point. 
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